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ABSTRACT 


Rihaczek  has  defined  the  dynamic  bandwidth  of  a signal 
as  the  square  root  of  the  magnitude  of  the  derivative  of  the 
instantaneous  frequency.  This  tutorial  report  reviews  the 
work  of  Rihaczek  and  Baghdady  on  the  practical  applications 
of  this  quantity.  It  is  shown  that  the  dynamic  bandwidth  is 
useful  as  a measure  of  1)  the  frequency  spread  of  the  signal 
energy  about  the  instantaneous  frequency,  2)  the  time  sta- 
bility of  the  instantaneous  frequency,  and  3)  the  minimum 
bandwidth  of  a filter  that  can  pass  the  signal  without  signif- 
icant phase  distortion. 

V 


_ 


_ 


GENERAL  $$  ELECTRIC 


r 


Page 


1.  INTRODUCTION 1 

2.  SIGNAL  ENERGY  DISTRIBUTION 2 

3.  FILTER  BANDWIDTH 6 

4.  EXAMPLES 8 

REFERENCES 11 

APPENDIX  A 13 

APPENDIX  B 15 


GENERAL  ■$&  ELECTRIC 


dynamic:  bandwidth 


T.G.  Kincaid 


1.  INTRODUCTION 

The  purpose  of  this  report  is  to  bring  attention  to  the  measure  of  "dynamic 
bandwicth"  defined  by  Rihaczekd)  and  to  urge  its  use  as  an  aid  in  signal  model- 
ing and  system  design. 

Consider  a complex  signal. 


x(t)  = a(t)  exp[j9(t)]. 


(1.  1) 


where  the  amplitude  a(t)  and  the  phase  9 (t)  are  real.  The  well-known  defini- 
tion of  the  instantaneous  frequency  of  x(t)  is 


- 2 


1 d Q (t ) 


(1.  2) 


where  f^(t)  is  in  hertz  when  t is  in  seconds.  Less  well  known  is  the  dynamic 
bandwidth  of  x(t)  defined  by 

df;(t)  i ll2 

(1.3) 


where  b(t)  is  in  Hertz  when  t is  in  seconds.  Note  that  the  dynamic  bandwidth 
is  related  to  the  phase  by 


. 1 d29  1/2 

b(t)  - 2rr  d? 


(1.  4) 


The  rationale  for  the  definition  of  dynamic  bandwidth  is  contained  in  the 
work  of  Rihaczek,  ^ and  the  definition  is  further  justified  by  the  earlier  work 
of  Baghdady.  '2,  3)  results  obtained  by  Rihaczek  (discussed  in  Section  2) 

show  that  the  dynamic  bandwidth  has  meaning  at  each  time  point  in  the  sense 
of  the  signal  energy  being  concentrated  in  a range  of  frequencies  about  the 
instantaneous  frequency.  He  also  shows  that  the  reciprocal  of  the  dynamic 
bandwidth  is  a measure  of  the  time  over  which  the  instantaneous  frequency  is 
essentially  constant.  The  work  of  Baghdady  (discussed  in  Section  3)  relates 
the  dynamic  bandwidth  to  the  narrowest  filter  bandwidth  that  can  pass  the  sig- 
nal without  significant  phase  distortion. 

A word  about  what  the  dynamic  bandwidth  does  not  measure  is  advisable. 
First,  since  the  dynamic  bandwidth  takes  only  frequency  fluctuations  into  ac- 
count, it  does  not  include  any  amplitude  modulation  contributions  to  the  signal 
bandwidth.  It  is,  therefore,  most  useful  for  signals  with  slow  amplitude  mod- 
ulation (compared  to  the  frequency  modulation).  Second,  the  dynamic 
bandwidth  is  not  a measure  of  the  overall  signal  bandwidth.  It  only  measures 
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the  frequency  spread  at  a particular  instant.  Large  slow  swings  in  instanta- 
neous frequency  can  result  in  an  overall  signal  bandwidth  much  larger  than 
the  maximum  dynamic  bandwidth.  Third,  when  the  dynamic  bandwidth  is  used 
to  specify  the  bandwidth  of  a low  phase  distortion  filter,  the  low  phase  distor- 
tion occurs  at  all  frequencies,  but  significant  amplitude  distortion  can  occur 
unless  the  signal  is  confined  to  the  filter  passband. 

In  summary,  the  results  described  in  this  report  show  that  the  dynamic 
bandwidth  gives  the  following  information  at  each  time  point: 

1.  The  frequency  spread  of  the  signal  energy  about  the  instantaneous 
frequency. 

2.  The  time  interval  over  which  the  instantaneous  frequency  is  stable. 

3.  The  minimum  bandwidth  of  a filter  which  can  pass  the  signal  without 
significant  phase  distortion. 

2.  SIGNAL  ENERGY  DISTRIBUTION 

Many  authors  [see  references  in  Rihaczek^]  have  described  signal  energy 
as  simultaneously  distributed  in  time  and  frequency.  The  work  of  Rihaczek 
collects  the  ideas  of  the  previous  authors  into  a unified  theory.  For  a com- 
plex signal  x(t),  with  Fourier  transform  X(f),  he  summarizes  the  previous 
descriptions  by  defining  the  complex  energy  density  in  time  and  frequency  as 


e(t,  f)  = x(t)  X(f)  exp[-j2rrft] 


(2.  1) 


He  justifies  this  definition  by  showing  that  1/2  the  integral  of  e(t,  f)  within  a 
time  interval  T and  frequency  band  B does  indeed  give  the  signal  energy  within 
this  time -frequency  cell. 

Rihaczek' s investigation  of  the  interesting  properties  of  the  complex  en- 
ergy density  function  is  recommended  reading,  but  contains  more  information 
than  is  necessary  here.  For  the  purposes  of  this  discussion,  only  his  use  of 
the  function  to  examine  the  distribution  of  energy  in  time  and  frequency  is 
important.  To  begin  this  investigation,  he  first  writes  the  signal  and  its 
Fourier  transform  in  terms  of  their  amplitude  and  phase;  i.  e. , 


x(t)  = | x(t)|  exp[j0(t)] 
X(f)  = |X(f)|exp[  j0  (f)] 


(2.2) 


(2.  3) 


If  these  are  substituted  into  Eq.  (2.  1),  then  the  signal  energy  E is  shown  in 
(1)  to  be  given  by  the  integral 

E = 1/2 j"J dtdf  | x<t)HX(f)  ] exp[j0(t)-j0(f)-2rrft]  (2.4 
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The  time  and  frequency  integrals  above  can  each  be  evaluated  by  the  method 
of  .stationary  phase  (see  Appendix  A).  This  method  uses  the  fact  that  the  sig- 
nificant contributions  to  the  integrals  occur  at  the  points  where  the  phase  of 
the  integral  is  stationary,  i.  e.,  the  phase  function  has  zero  slope.  The  method 
requires  that  the  second  derivative  of  the  phase  is  not  zero  at  these  points,  and 
that  the  amplitude  variation  is  slow  compared  to  the  phase  variation  in  the  re- 
gion about  the  points.  The  stationary  points  can  be  found  by  setting  the  time 
and  frequency  phase  derivatives  equal  to  zero,  and  solving  for  the  coordinates 
of  the  stationary  points. 


Setting  the  time  derivative  of  the  phase  of  e(t,  f)  equal  to  zero  gives  the 
following  equation  for  the  times  at  which  the  phase  is  stationary: 


f - 


1 dg(t) 

2tt  dt 


(-fi(O) 


(2.  5) 


This  equation  is  just  the  definition  of  the  instantaneous  frequency  of  x(t).  For 
each  f,  the  values  of  t which  solve  this  equation  are  the  stationary  time  points 
corresponding  to  frequency  f.  For  any  physical  signal,  x(t)  has  an  instanta- 
neous frequency  over  its  entire  time  extent;  thus  every  time  point  on  the  ex- 
tent of  x(t)  is  a stationary  time  point  for  some  f. 


Setting  the  frequency  derivative  of  the  phase  of  e(t,  f)  equal  to  zero  gives 
the  following  equation  for  the  frequencies  at  which  the  phase  is  stationary. 


t 


J_  d»tf) 

2tt  df 


(=  tg(f)) 


(2.6) 


This  equation  is  the  definition  of  the  familiar  group  delay  tp  of  x(t),  the  dual 
of  the  instantaneous  frequency.  At  each  t,  the  values  of  f which  solve  this 
equation  are  the  stationary  frequency  points  corresponding  to  that  particular 
value  of  the  time  t.  For  any  physical  signal,  x(t)  has  a group  delay  for  all 
frequencies  in  its  band.  Therefore,  every  frequency  in  the  band  occupied  by 
x(t)  is  a stationary  point  for  some  t.  (Note  that  for  theoretical  signals  with 
spectral  line  components,  the  group  delay  is  not  defined  at  the  line  frequencies.  ) 

It  has  been  argued  that  almost  every  value  of  t and  f within  the  time -fre- 
quency extent  of  x(t)  are  stationary  points  for  their  respective  integrals. 
However,  a significant  contribution  to  the  double  integral,  and  hence  to  the 
signal  energy,  occurs  only  in  the  region  of  the  t,  f plane  where  the  stationary 
points  are  simultaneous  solutions  to  Eqs.  (2.  5)  and  (2.  6).  Each  of  these 
equations  describes  curves  in  the  time  frequency  plane.  Therefore,  their 
simultaneous  solutions  are  at  the  intersections  of  these  two  curves.  It  is 
shown  in  Appendix  B that  these  two  curves  intersect  all  along  their  length,  i.  e. , 
they  are  the  same  curve.  Mathematically,  this  means  that  f^( ) and  tg( ) are 
functional  inverses.  That  is,  * 

fi(tg(f>)  = f (2.  7) 

Therefore,  either  Eq.  (2.  5)  or  (2.  6)  is  the  locus  of  all  the  points  in  the  t,  f 
plane  at  which  e(t,  f)  is  simultaneously  stationary  in  both  t and  f.  It  follows 
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that  at  each  point  in  time  the  signal  energy  is  concentrated  in  an  inteval  about 
the  instantaneous  frequency.  The  width  of  this  interval  in  the  frequency  direc- 
tion is  the  band  of  frequencies  over  which  the  phase  of  e(t,  f)  is  nearly  a con- 
stant as  a function  of  frequency.  Similarly,  the  width  in  the  time  direction  is 
the  duration  of  time  over  which  the  phase  of  e(t,  f)  is  nearly  constant  as  a 
function  of  time.  It  is  shown  below  that  these  widths  are,  respectively,  the 
dynamic  bandwidth  and  a quantity  called  the  signal  relaxation  time. 


Rihaczek  determines  the  width  of  these  intervals  by  noting  that  the  re- 
quirement for  the  phase  of  e(t,  f)  to  be  nearly  constant  on  an  interval  is  that 
9(t)  and  0(f)  are  nearly  linear  on  the  interval.  He  defines  these  intervals  to 
be  those  over  which  the  quadratic  terms  in  the  Taylor  expansions  of  0 (t ) and 
0(f)  do  not  deviate  from  linearity  by  more  than  tt/4.  He  uses  this  criterion  to 
determine  the  time  interval  width  by  noting  that  the  magnitude  of  the  quadratic 
term  in  the  Taylor  expansion  of  8(t)  about  the  stationary  point  t0  is 


1/2  (t-t0)2 


d20 


dt2 


By  using  this  term  as  an  approximation  to  the  deviation  from  linearity,  the 
interval  can  be  determined  as  the  time  between  the  two  points  where  its  mag- 
nitude is  tt/4.  Computing  this  number,  called  the  signal  relaxation  time  r, 
gives 

(2.8) 


1 

i/2 

1 

1 d20(t) 

df^t) 

2 it  dt2 

dt 

when  the  derivatives  are  evaluated  at  time  t0.  Using  the  same  criterion,  the 
frequency  interval  about  the  stationary  point  f0  over  which  the  phase  is  nearly 
constant  can  be  shown  to  be 


1 

1/2 

1 

1 d20  (f) 

dtg(f) 

2tt  df2 

df 

(2.  9) 


when  the  derivatives  are  evaluated  at  fo.  The  use  of  b is  deliberate,  since 
this  is  the  dynamic  bandwidth  defined  in  Eq.  (1.  3).  This  is  proved  from  the 
fact  that  f j and  tg  are  inverse  functions,  and  hence  their  derivatives  (i.  e. , 
their  slopes)  are  reciprocals  at  the  same  point  on  their  trajectory.  That  is. 


idf  i (t) 


dt 


dtg(f) 


df 


when  evaluated  at  the  same  point  on  the  projectory  of  fj(t)  and  tg(f). 
relationship  is  pictured  in  Figure  2.  1. 

Substitution  of  Eq.  (2.  10)  into  Eq.  (2.  9)  gives 

1/2 

.1  ~A  I ' I * w OVW  I • 


(2.  10) 


This 


dfi<t)|i/2 

1 d29(t) 

dt  \ 

2tt  dt2 

(2. 11) 
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trajectory  of 
fj(t)  and  tg(f) 


Slope  = 


df^t) 


dtg(  f) 

df 


Figure  2.1.  The  slope  of  the  instantaneous  frequency  and  group  delay 
trajectory  at  the  point  t0,  fo  in  the  time -frequency  plane. 

which  is  the  definition  of  dynamic  bandwidth  given  in  Eq.  (1.  3).  A comparison 
of  Eqs.  (2.  8)  and  (2.  11)  shows  that  b and  r are  reciprocals  when  computed  at 
the  same  point  on  the  fj(t)  and  tg(f)  trajectory.  This  leads  to  the  time- 
frequency  uncertainty  relation  derived  by  Rihaczek, 


br  = 1. 


(2.  12) 


This  result  shows  that  at  each  instant  of  time  the  signal  energy  is  concentrated 
in  a unit  area  about  the  instantaneous  frequency  trajectory  in  the  time- 
frequency  plane.  This  area  changes  in  shape  as  time  progresses.  Slow  in- 
stantaneous frequency  changes  spread  the  energy  more  in  time,  while  rapid 
changes  spread  it  more  in  frequency.  These  concepts  are  pictured  in 
Figure  2.2. 
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Figure  2.2.  Signal  energy  concentration  about  the  instantaneous  frequency 
trajectory  at  two  points  in  the  time  frequency  plane. 


3.  FILTER  BANDWIDTH 

Baghdady^»  3)  has  determined  the  minimum  filter  bandwidth  consistent 
with  low  distortion  filtering  of  a frequency  (i.  e.,  phase)  modulated  signal. 

He  determines  this  by  considering  a system  in  which  a pure  phase  modulated 
signal,  x(t)  = exp[j9(t)J,  is  the  input  to  a filter  with  radian  frequency  response 
H(ju).  The  filter  output  y(t)  has  "low  distortion"  if  at  each  time  instant  y(t)  is 
approximately  the  response  of  the  filter  to  a single  frequency  input  at  the  in- 
stantaneous frequency  of  x(t).  That  is,  the  filter  has  "low  distortion"  if 


y(t)  =“  H(0'(t))  exp[j9(t)] 


(3.  1) 


where  the  prime  denotes  derivative.  Baghdady  shows  that  the  actual  response 
is  given  by 
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y(t)  II(o'(t))  ll+e(t)]  explj9(t)J 


(3.2) 


where  c(t)  is  a relative  error  term  which  satisfies  the  condition 

I e (t) ' ^ 1/2 1 9"  (t) ' max  I^TTT 

1 max  | u(fjj)  max 


(3.  3) 


For  lov/  distortion,  the  right  side  of  the  above  equation  must  be  « 1;  i.  e., 
the  relative  error  must  be  small. 

The  ratio  I H"  (uu)/H(x)  | is  called  the  sluggishness  ratio.  It  is  a measure 
of  the  ability  of  the  energy  storage  elements  in  the  filter  to  respond  to  changes 
in  the  instantaneous  frequency  of  the  signal.  Baghdady  has  calculated  this 
ratio  for  Butterworth  filters  up  to  sixth  order,  (2,  3)  and  for  cascaded  first  - 
and  second -order  Butterworth  filters  up  to  six  sections.  (3)  In  all  cases,  the 
maximum  value  of  the  sluggishness  ratio  is  inversely  proportional  to  the 
square  of  the  filter  half  power  bandwidth.  This  relation  appears  to  be  gener- 
ally true,  and  the  fact  that  it  is  true  for  many  filters  of  interest  allows  the 
maximum  value  of  the  ratio  for  these  filters  to  be  written 


III"  (ju) 
H (yu) 


k 

(2ttBF 


(3.  4) 


The  bandwidth  B is  in  Hertz,  and  the  proportionately  constant  k is  called  the 
stiffness  index  of  the  filter. 

Baghdady  gives  the  value  of  k for  the  filters  mentioned  above.  (3)  For  a 
first-order  Butterworth  filter  (i.  e.,  a single  tuned  resonant  circuit)  k = 8. 

The  index  increases  for  higher  order  Butterworth  filters.  A simple  calcula- 
tion shows  that  for  a sine  function  filter  k = tt2  if  the  bandwidth  is  defined  as 
the  distance  from  the  peak  of  the  spectrum  to  the  first  zero.  An  interesting 
case  is  the  filter  with  an  exponentially  decaying  frequency  response  exp[-|f  |/B], 
for  which  k = 1 if  the  bandwidth  is  defined  as  the  distance  from  the  peak  to  the 
frequency  where  the  amplitude  is  down  by  l/e. 

Substituting  Eq.  (3.  4)  and  the  definition  of  instantaneous  frequency 
[Eq.  (1.  4)]  into  Eq.  (3.  3)  gives  an  alternate  expression  for  the  relative  error 
bound. 


. . k b2max 

!e<t>  5 TZ  “52— 


(3.5) 


For  low  distortion  reproduction  of  the  phase  modulation  at  the  filter  output, 
the  right  side  of  the  above  inequality  should  be  on  the  order  of  1/10.  For 
convenience,  make  it  l/4n.  Then  the  minimum  bandwidth  of  the  low  distortion 
filter  is  given  by 


bT 


(3.  6) 


This  equation  shows  that  the  filter  bandwidth  is  proportional  to  the  signal  dy- 


namic bandwidth,  as  it  should  be  for  a meaningful  definition  of  dynamic 
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bandwidth.  The  constant  of  proportionality  depends  upon  the  type  of  filter 
through  its  stiffness  index  k.  For  a first-order  Butterworth  filter,  k 8 and 
B - - bniax*  *’or  a sine  filter,  k and  B - bmax.  For  the  exponential 
filter  mentioned  above,  k 1 and  B bmax- 

A word  about  the  practical  uses  of  Eq.  (3.  5)  to  specify  filter  bandwidth. 
This  design  equation  was  derived  to  meet  a low  phase  distortion  criterion. 

The  low  phase  distortion  will  hold  at  all  frequencies,  even  outside  the  filter 
pass  band.  However,  the  amplitude  of  the  signal,  in  general,  will  be  distorted 
outside  the  pass  band.  If  Eq.  (3.5)  is  used  to  design  a bank  of  filters--  say, 
for  spectrum  analysis  - -the  criterion  given  by  Eq.  (3.  5)  is  ideal  for  determin- 
ing the  minimum  bandwidth,  since  the  signal  will  always  be  in  one  of  the  filters. 
However,  for  a tuning  or  noise  rejection  filter,  care  must  be  taken  that  the 
signal  remains  in  the  pass  band  to  prevent  amplitude  distortion. 


4.  EXAMPLES 


4.  1 Single  Frequency 


x(t)  = exp[j2rrf0t] 

Phase:  9(t)  = 2TTf0t 

Instantaneous  Frequency:  f£(t)  = f0 

Dynamic  Bandwidth:  b (t ) = 0 

The  zero  dynamic  bandwidth  is  consistent  with  the  notion  of  zero  bandwidth 
for  a single  frequency  signal.  Note  that  the  group  delay  and  relaxation  time 
are  not  defined  for  this  signal.  The  instantaneous  frequency  trajectory  of  this 
signal  is  just  a straight  line  parallel  to  the  time  axis,  as  shown  in  Figure  4.  1. 


Figure  4.  1.  The  instantaneous  frequency  trajectory  for  a single 
frequency  signal. 
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4.  2 r M Slide 

x(t)  expfj^at2] 

Phase:  6(0  -trat2 

Instantaneous  Frequency:  fj[(t)  at 

Dynamic  Bandwidth:  b(t)  - | a.  | 1 / 2 

The  instantaneous  frequency  of  the  FM  slide  has  a linear  time -frequency 
trajectory  and  a constant  dynamic  bandwidth.  The  group  delay  is  f/a  and  the 
relaxation  time  is  | a j _1/2.  The  instantaneous  frequency  trajectory  is  shown 
in  Figure  4.2,  along  with  the  frequency  distribution  of  energy  given  by  the  dy- 
namic bandwidth. 


Figure  4.2.  The  frequency  distribution  of  signal  energy  about  the  instanta- 
neous frequency  trajectory  of  an  FM  slide. 
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4.  3 Sinusoidal  Phase  Modulation 


Phase: 

Instantaneous  Frequency: 
Dynamic  Bandwidth: 


x(t)  = exp[jgsin2rrf0t] 

9(t)  = 0sin2rrfot 

l'i(t)  = pf0cos2rrf0t 

b(t)  f0 1 2tt  3sin2rrf0t  | 1 /5 


When  the  phase  modulation  is  sine  wave  the  instantaneous  frequency  is  a 
cosine  of  the  same  frequency,  and  the  dynamic  bandwidth  oscillates  at  twice 
the  modulation  frequency.  Since  the  signal  x(t)  is  periodic,  it  has  a line  spec- 
trum, and  the  group  delay  and  relaxation  time  are  not  defined. 

The  instantaneous  frequency  trajectory  is  shown  in  Figure  4.  3,  along  with 
the  frequency  distribution  of  energy  given  by  the  dynamic  bandwidth.  The 
maximum  dynamic  bandwidth  is  /2rrp  fo,  and  occurs  when  the  frequency  modu- 
lation (i.  e.,  instantaneous  frequency)  goes  through  zero.  This  maximum 
value  determines  the  minimum  bandwidth  of  a filter  required  to  pass  x(t)  with- 
out phase  distortion.  Note  that  there  can  be  amplitude  distortion  if  the  am- 
plitude of  the  instantaneous  frequency  pf0  is  large  compared  to  the  dynamic 
bandwidth. 


Figure  4.  3.  The  frequency  distribution  of  signal  energy  about  the  instanta 
neous  frequency  trajectory  of  a sinusoidally  phase  modulated 
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APPENDIX  A 


Tnis  Appendix  is  a summary  of  the  method  of  stationary  phase  as  given  by 


Papoulis.^  5 

Consider  the  integral 


uu 

I J a(t)  exp[ja(t)] 


(Al) 


If  a(t)  varies  rapidly  compared  to  a (t),  then  I can  be  evaluated  by  the  method 
of  stationary  phase.  Assume  a!(t)  - 0 at  only  one  point  to,  called  the  stationary 
point  of  alt ).  Then 


I » a 


(to)/i3^exp[ja(to)± 


(A2 ) 


w 


here  the  sign  in  front  of  tt/4  is  the  same  as  the  sign  of  a"  (t0). 


Suppose  that  a(t)  has  n stationary  points,  ti,  t2,  — , tn.  Then  the  method 
of  stationary  phase  gives 

<A3) 

where  the  sign  of  r/4  is  determined  as  above. 
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APPENDIX  B 


Consider  the  signal 


x(t)  = |x(t)|  exp[j9(t)j 


(Bl) 


on  the  interval  T.  The  instantaneous  frequency  of  x(t)  is  the  function 


fi(t)  = ^ 


d G (t ) 
dt 


(B2) 


Assume  the  interval  T is  chosen  so  that  no  value  of  the  instantaneous  frequency 
occurs  more  than  once.  Mathematically,  this  means  that 


dfi(t)  d2e(t)  , 

dt  “ dt*  f 


(B3) 


on  the  interval  T. 

The  Fourier  transform  of  x(t)  on  T is  given  by 

X(f)  = J dt|x(t)|[exp  j9 (t)  - j2nft] 

T 

= lX(f)|  exp[j0(f)]  (B4) 

The  group  delay  of  x(t)  on  T is  the  function 

(B5) 

The  instantaneous  frequency  is  a function  of  time,  and  the  group  delay  is  a 
function  of  frequency.  Each  of  these  functions  describes  a trajectory  in  the 
time -frequency  plane.  It  is  shown  below  that  these  two  trajectories  are  the 
same;  i.  e. , f^(  ) and  tg(  ) are  inverse  functions. 

The  proof  of  this  is  as  follows.  For  any  value  of  f,  the  stationary  point 
of  the  signal  phase  is  the  time  point  ts  which  solves  the  equation 

f = JLlilll  (B6) 

2tt  dt 

Choose  f so  that  tg  is  not  near  the  edge  of  the  interval.  Then,  by  the  principal 
of  stationary  phase,  the  Fourier  transform  X(f)  is  closely  approximated  by 

X(f)=“  ]x(t)|  exp[j0(t)  - j2nft  ± Tr/4]|ts  (B7) 

where  the  sign  of  rr/4  is  chosen  the  same  as  the  sign  of  9 " (ts ).  Thus,  the 
phase  of  X(f)  is 

0(f)  = B(ts)  - 2ntsf  ± rr/4  (B8) 
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Note  that  t is  a function  of  f.  By  definition 

t (f) , . jl 

g'  2tt  df 


1 d9(ts)  . * , dts  . 

2rr  df  S + df  f 


1 dts 


— + t + — — f 

2tt  dtg  df  ls  df 


However,  since  tg  is  the  solution  to  Eq.  (B2)  for  f^tg)  f;  therefore. 


d0(ts) 


- 2 TT  f 


(BIO) 


Substitution  of  (B9)  into  (B8)  gives 


tg(f)  = ts(f) 


(Bll) 


That  is,  the  stationary  time  point  for  any  f is  the  group  delay  time  for  that  f. 
Therefore,  tg(f)  is  a solution  to  Eq.  (A2),  and 


W»  = 2 


(B12) 


which  was  to  be  proved. 
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